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Abstract - This article describes and analyses the existing 
possibilities for using Standard Statistical Methods and 
Artificial Intelligence Methods for a short-term forecast and 
simulation of demand in the field of Asymmetric Digital 
Subscriber Line Internet Connection in the Czech Republic. 
The most widespread methods are based on Time Series 
Analysis. Nowadays, approaches based on Artificial 
Intelligence Methods, including Neural Networks, are 
booming. Separate approaches will be used in the study of 
Demand Modelling in the field of Asymmetric Digital 
Subscriber Line, and the results of these models will be 
compared with actual guaranteed values. Then we will examine 
the quality of Neural Network models. The another part of 
study will be focused on improving the quality of Neural 
Network models with the use of indicator Gross Domestic 
Product. 

Index Terms - Demand, Telecommunications, Standard 
Statistical Methods, Neural Network, Asymmetric Digital 
Subscriber Line, Gross Domestic Product 

I. INTRODUCTION 

Demand can be defined as the relation between price and 
the quantity of goods that buyers are willing to purchase. 
This correlation is displayed in relation to the global market 
by the sold product quantity at one time-point. If we focus 
on the telecommunication services sector, we can note the 
development of the sale of cell phones and internet 
connection (Asymmetric Digital Subscriber Line, Wireless 
connection, etc.). This article will be focused on Demand 
Modelling of the Asymmetric Digital Subscriber Line internet 
connection (ADSL) in the Czech Republic. Considering the 
fact that in the literature listed in the paper, an effective 
procedure of demand modelling of the Asymmetric Digital 
Subscriber Line Internet Connection was not described, this 
research project conceives the demand as time series under 
which the demand model can be designed and trends 
predicted. The demand model was formed on the basis of the 
data of statistical survey on the territory of the Czech Republic 
(e.g. number of households in the Czech Republic who have 
the technology of ADSL) which were published by the Czech 
Statistical Office in the period of 2006 - 20 1 1 . The theoretical 
description of the methods for construction of the demand 
model is listed below. 

H. CONSTRUCTION OF THE DEMAND MODEL 

The demand model can be constructed using standard 
statistical methods including Decomposition Time Series, and 
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the Box-Jenkins Metodology, or by applying Artificial 
Intelligence methods, including for example Neural Networks 
[8]. 

A. Decomposition Time Series 

The series {y, t= 1, T} is gradually decomposed to 
several components: trend, circular component, seasonal 
component and residual component (unsystematic 
component). This method is based on work with time series 
systematic components. Features of time series behaviour 
can be better observed in separate components than in the 
undecomposed original time series [8]. In this research study 
from the field of standard statistical methods, exponential 
smoothing will be used for demand modelling of the 
Asymmetric Digital Subscriber Line internet connection. 

/. Exponential Smoothing 

The above defined time series will be written as { y , 
t = 1, T}. Simple Exponential Smoothing is described in the 
recurrent form 

y t =ay t +(l-a)y M , 

y is the Exponential Average in time t, y _ x is the 

Exponential Average in time t- 1 , value a is the Smoothing 

Constant from the interval OC G< 0;1 > [8]. The Exponential 

Average can be expressed on the basis of the recurrent form 
as [8]: 

9, =ay,+(\-a)y,-i =«y, + (l-«)[«y,_, + (l-a)y,_ 2 ] = 

= ay, +a(l-a)y M +(l-a) 2 [ay t _ 2 (l-a)y t _ 3 ] = ...= 

= ay, +a(l-a)j,_, +{l-a) 1 ay l _ 2 +... + a(l-a)' y t _, +.... 

.. + (1 - a)' y = a£(l - a)' y,_, + (1 - a)' % 

(=0 

//. Brown 's Simple Exponential Smoothing 

The time series y t is constructed with a stationary process 

in the form y t = /3 Q + E t , ft is the mean value of the 

process, and S t are random values with the features of white 

noise. After applying Exponential Smoothing to the Time 
Series we obtain the relation [8] : 

y, =af)l-a)';y M =a%l-aj(fi ) +£,_,)=#, +a^}l-a)' £,_, 

i=0 i=0 i=0 
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because 

a£(l-a) = l 

applies to the mean value and 

dispersion: 

E(y c ) = E(y t ) = /fc f = = — ^ 

2 — a, 2 — oc ' 
Brown's Linear (double) Exponential Smoothing, alternatively, 
we can make a double application of the Simple Exponential 
Smoothing method on the time series y expression in the 
form [8]: 

B. Box-Jenkins Methodology 

Unlike classical decompositional methods, which deal 
with systematic time series components (trend, circular and 
seasonal), the Box-Jenkins methodology deals with a 
residual (unsystematic) component. The method involves 
searching for relations of individual observations. By this 
method we are able to describe time series which are not 
manageable by standard methods [8]. The time series is 
perceived as a realization of the stochastic process which is 
defined as a series of accidental quantities arranged in time 

[X (j,(),seS,(£T}, where S is a selective space and T is an 
index series. For each s S the realization of the stochastic 
process is defined on the index series T [8]. For the Box- 
Jenkins methodology, the following special accesses are 
significant: autoregressive process AR, moving average 
process MA, and combined process ARMA. These 
processes result from the linear process by resetting all 
parameters till the final number. The parameters are chosen 
in such a way that the stationarity and invertibility of the 
processes will be ensured. A special non-stationary ARIMA 
model also exists in the Box-Jenkins metodology [8]. 

/. ARIMA Processes 

Some integrated processes maybe arranged by means of 
differentiation to stationary and are expressed in the form of 
the stationary and invertible ARMA(p,</) model. The original 
integrated process in the form [8]: 

<P p (B)(l-B) d y, =6 q {B)s, 

is called the autoregressive integrated process of sliding 
averages of the order p, d, q. This is called ARIMA (p, d, 
q). Models with d = 1,2 are usually used [8]. 

C. Artificial Intelligence Methods - Neural Networks 

The benefit of the Artificial Neural Network lies in its 
ability to implement complex non-linear functions. Neural 
systems co-execute a large number of operations and work 
without an algorithm [8]. Their activity is based upon the 
learning process, when the neural network gradually 
conforms to calculation. In the course of a learning phase 
we do not have to be occupied by the problem of the right 
selection function, because the neural network is able to 
make do only with practise examples [8]. This is the main 
difference in comparison with a traditional approach (e.g. 
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comparison of traditional non-linear models and a back- 
propagation network). The neurone function scheme can be 
demonstrated graphically as follows [8]: 



n*miin*ar 

Tyn? lion 

/ 




Figure 1. Model of a neurone [8] 

N 

y = S y^w^ +© 

_ w 

where y is the output (neurone activity), S demonstrates a 
transmission function (jump, linear, non-linear), X is neurone 

input (inputs are in total N), W, represents a level of synoptic 

weight, describes a trigger level. 

I. The Principle of Back-Propagation Network Learning 
Let us consider a neurone network with L - layers 1=1, 
. . ..,L and as the output i-th neurone in the 1-th layer we use 

the indication V- . V° means X t , i.e. the i-th output. The 

indication w\j expresses a connection from Vf 1 to . V- . 

An algorithm can then be inscribed after separate stages as 
follows [8] : 

1. Subranged random numbers based - Weight 
Initialization 

2. Insertion of — p into the network input (layer 1 = 0), i.e. 

3. Network signal propagation: 

v! = g(h!) = 8 (^w l ij v;- 1 ) 

k 

4. Calculation of delta for the output layer: 

5. Calculation of delta for previous layers by Error Back- 
Propagation 

S? =E'(hr^ l ,8' i 



6. Weight change according to formula: 

/ „ e Zt t I— 1 new old 



H =r }8;v;- 



new 

W- = W- 

B B 



+ Aw,,. 



5+ 



7. If all samples have been submitted to the network 
we continue in phase no. 8, otherwise we go back to 
phase no. 2. 
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8. If the network error compared to the selected criterion 
value was minor or the maximum number of steps was 
exhausted, then the learning process can be completed, else 
phase no. 2. 

D. Choice of a Relevant Model 

An appropriate model can be determined on the basis of 

a) the graph of the time series, or from its absolute or relative 
characteristics [8], 

b) interpolative criteria (a decisive deviation of residues, 
coefficient of determination, coefficient of autocorrelation of 
residues, average characteristics of residues ) [8], 

c) extrapolative criteria (average characteristics of "ex post" 
forecast mistakes, graph forecasts) [8]. 

E. Average Characteristics of Residues 

The average square mistake - dispersion [8]: 

MSE = ±±(y t -y t f=±-±a? 



n — 



The root mean square mistake [8]: 



RMSE= ±±(y-%f = ±±a? 
\n l=l \n t=i 

The average absolute mistake [8]: 



MAE 



n Z 



5>. 



i - 

n 



The lower the values of the specified characteristics, the 
better the chosen model is [8]. 

III. DEMAND MODELLING OF ASYMMETRIC 
DIGITAL SUBSCRIBER LINE INTERNET CONNECTION 

The above mentioned methods will be used for demand 
modelling of the Asymmetric Digital Subscriber Line Internet 
Connection (ADSL) in the Czech Republic. The ADSL Demand 
Model (i.e. the number of households which are using a 
Asymmetric Digital Subscriber Line Internet Connection) is 
based on the values from the Czech Statistical Office (between 
2006 and 201 1) and estimate made by an expert. The results 
for the ADSL demand model based upon Box-Jenkins 
methodology, Brown's linear exponential smoothing and 
Neural Networks are graphically displayed below. 

A. ARIMA Processes 

The figure 2 describes the comparison of the real values 
(i.e. the number of households which are using an 
Asymmetric Digital Subscriber Line Internet Connection), and 
the values from a demand model based on AREVIA (1,1,0). 



Demand Modelling of ADSL 
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Figure 2. ARIMA Demand Modelling of ADSL 

The figure 3 shows the course of residual values (i.e. the 
difference between the number of Households which are 
using the Asymmetric Digital Subscriber Line Internet 
Connection, and the value predicted by demand model based 
on ARIMA( 1,1,0)). 
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Figure 3. Residual values for ARIMA Demand Modelling of ADSL 

In the figure 4, we can observe the course of residual 
autocorrelations for a demand model based on ARIMA( 1 , 1 ,0). 
In this case, the residual autocorrelation was not recorded. 



Estimate Autocorrelations for ADSL Demand Model 
ARIMA(1,1,0) 
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Figure 4. Residual Autocorrelations for ADSL 

B. Brown 's linear exponential smoothing 

The figure 5 shows the comparison of the real values (i.e. 
the number of households which are using an Asymmetric 
Digital Subscriber Line Internet Connection), and the values 
from a demand model based on Brown's linear exponential 
smoothing. 
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Demand Modelling of ADSL 
Brown's exp. smoothing with alpha = 0,7794 
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Figure 5. Brown's linear exponential smoothing of ADSL 

In the figure 6, we can observe the course of residual 
autocorrelations for a demand model based on Brown's lin- 
ear exponential smoothing. In this case, the residual 
autocorrelation was not recorded. 



Estimate Autocorrelations for ADSL Demand Model 
Brown's exp, smoothing with alpha = 0,7794 
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Figure 6. Residual Autocorrelation for ADSL 
C. Back-Propagation Neural Network 

The figure 7 shows the comparison of the real values (i.e. 
the number of households which are using the Asymmetric 
Digital Subscriber Line Internet Connection), and the values 
from a demand model based on Neural Network. 



Demand Modelling of ADSL 
Neural Network 
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Figure 7. Neural Network Demand Modelling 

The course of residual values (i.e. the difference between the 
number of Households which are using the Asymmetric Digi- 
tal Subscriber Line Internet Connection), and the value pre- 
dicted by demand model based on Neural Network is shown 
in figure 8. 
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Residual values for ADSL 
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Figure 8. Residual values for ADSL 

The figure 9 shows the course of residual autocorrelations 
for a demand model based on Neural Network. In this case, 
the residual autocorrelation was recorded. 



Estimate Autocorrelations for ADSL Demand Model 
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Figure 9. Residual Autocorrelation for Neural Network 

Within the research project on demand modelling in the field 
of Asymmetric Digital Subscriber Line Internet Connection 
in the Czech Republic we obtained the following values for 
the average characteristics of residues, as shown in Table I. 

Table I. Average Characteristics Of Residues 





Neural 
Network 


ARJMA. 
(1,1,0) 


Brown's exponential 
snicothuiE: 


EMSE 


5.24 


7.92 


1014 


MAE 


3.S9 


5.50 


6.02 



On the basis of the results of RMSE and MAE average char- 
acteristics of residues, the demand model based upon Neural 
Networks can be considered as best, but a different result 
arises from an evaluation of prediction quality, because the 
best prediction value result was obtained from a demand 
model based on ARIMA (1,1,0). In the figure 10, we can 
observe an evaluation of prediction quality for the ADSL 
demand model based upon ARIMA( 1,1,0), Brown's linear 
exponential smoothing and Neural Network. With respect to 
the above mentioned facts, another part of the study was 
focused on improving the quality or accuracy of the ADSL 
demand model with the use of neuron networks. 
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Demand Modelling In Telecommunications 

Comparison of Prediction Quality for ADSL 
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Figure 10. Comparison of the ADSL internet connection used by 
households in the 2nd quarter of 2010, values published on February 
10th, 2011 by the Czech Statistical Office, and results predicted 
by demand models for ADSL. 

D. Back-Propagation Neural Network (including GDP) 

Within the framework of the progressive proposal for 
bringing more precision to the current approach, the 
Asymmetric Digital Subscriber Line Internet Connection 
demand model was extended by the effect of the indicator 
Gross Domestic Product (GDP) of the Czech Republic and a 
Back-Propagation neural network was subsequently applied 
to the model formed in this way, as shown in figure 1 1 . 




In the figure 12, we can observe the comparison of the real 
values (i.e. the number of households which are using an 
Asymmetric Digital Subscriber Line Internet Connection), and 
the values from a demand model based on Neural Network 
which was extended by the effect of the indicator GDP Czech 
Republic. 



Demand Modelling of ADSL 
Neural Network (including GDP) 
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Figure 12. Neural Network (including GDP) Demand Modelling 

The course of residual values (i.e. the difference between the 
number of Households which are using the Asymmetric Digi- 
tal Subscriber Line Internet Connection, and the value pre- 
dicted by demand model based on Neural Network which 
was extended by the effect of the indicator GDP Czech Re- 
public) is shown in figure 13. 
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Residual values for ADSL 
Neural Network [including GDP) 
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Figure 13. Residual values for ADSL 

The figure 14 describes the course of residual 
autocorrelations for a demand model based on Neural 
Network (including GDP). The residual autocorrelation was 
not recorded. 



Estimate Autocorrelations for ADSL Demand Model 
Neural Network (including GDP) 
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Figure 14. Residual Autocorrelation for Neural Network (GDP) 

In this case, we obtained the following values for the aver- 
age characteristics of residues, as shown in Table II. 

Table II. Average Characteristics Of Residues 





Neural 
Xetwork (GDP) 


ARJMA 

cu,o) 


Brown's expcnentid 
smoothing 


RMSE 


431 


7 92. 


10 14 


MAE 


3.09 


5.50 


6.02 



On the basis of the results of RMSE and MAE average char- 
acteristics of residues, the demand model based upon Neural 
Networks (including GDP) can be considered as best and 
same result arises from an evaluation of prediction quality 
because the best prediction value result was obtained from 
a demand model based on Neural Network (including GDP), 
as shown in figure 15. 



Demand Modelling in Telecommunications 
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Figure 15. Comparison of the ADSL internet connection used by 
households in the 2nd quarter of 2010, values published on July 8th, 
2011 by the Czech Statistical Office, and results predicted by demand 
models for Asymmetric Digital Subscriber Line Internet Connection. 
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CONCLUSION 

This article deals with the demand modelling of the 
Asymmetric Digital Subscriber Line Internet Connection in 
the Czech Republic. Considering the fact that in the literature 
listed in the paper, an effective procedure of demand 
modelling of the Asymmetric Digital Subscriber Line Internet 
Connection was not described, this research project 
conceives the demand as time series under which the demand 
model can be designed and trends predicted. The 
introductory chapters summarize in theory individual 
approaches of formation of demand modelling based on 
statistical analysis of time series and further current 
procedures based on Artificial Intelligence Methods including 
Neural Networks. Furthermore, the objective of article was to 
evaluate contemporary possibilities concerning application 
of standard statistical methods and approaches of artificial 
intelligence with focusing on neural networks in forecasting 
and modelling demand in the segment of the Asymmetric 
Digital Subscriber Line Internet Connection in the Czech 
Republic. Individual demand models which were drawn up 
on the basis of the above described methods were evaluated 
and the results of forecast following from these models were 
compared with real development. On the basis of the results, 
it is possible to make a conclusion that the demand models 
formed with the aid of standard statistical methods prove 
higher quality or accuracy of the short-term forecast than the 
models drawn up by the aid of a neural network, as shown in 
figure 10. With respect to the above mentioned facts, another 
part of the study was focused on improving the quality or 
accuracy of the demand models with the use of neuron 
networks. Within the framework of the progressive proposal 
for bringing more precision to the current approach, the 
demand models in the segment of Asymmetric Digital 
Subscriber Line Internet Connection were extended by the 
effect of the indicator GDP Czech Republic and a neural 
network was subsequently applied to the model formed in 
this way, as shown in figure 11. The values resulting from 
such demand models proved positive effect of the indicator 
GDP on the quality or accuracy of the short-term forecast, as 
shown in figure 15. 



Demand Modelling of ADSL 
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Figure 16. Quality of prediction 
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With respect to the above mentioned results, the figure 16 
offers the percentage accuracy of the values predicted by 
demand models for Asymmetric Digital Subscriber Line 
Internet Connection which were drawn up on the basis of 
the above described methods. Unlike the approaches which 
have been applied so far, the accuracy of the demand models 
drawn up on the basis of neural networks was reached by the 
progressive method resting in implementation of the indicator 
GDP, as shown in figure 15, 16. The progressive method of 
modelling the demand in the course of which the results of 
the research project were obtained may potentially contribute 
to economic parameters of telecommunication companies, 
because quality demand models which provide precise results 
in a short-term forecast, contribute to telecommunication 
companies from the point of view of planning and evaluation 
of investment strategy, and furthermore, they positively affect 
the whole financial results and economic indicators of a 
company. 
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